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Abstract: We derive the mixture of odd to even barium isotopes in the atmosphere of the metal- 
poor subgiant HD 140283 from the analysis of the Ba ll transition at 4554 A in a high-resolution high 
signal-to-noise spectrum of the star. The detailed shape of this spectral line depends on the relative 
contributions of odd and even isotopes via isotopic and hyperfine splitting. We measure the fractional 
abundance of odd Ba isotopes by modelling the formation of the Ba ll 4554 A line profile with the use of 
both a classical ID hydrostatic and a 3D hydrodynamical model atmosphere of HD140283. We interpret 
the results in terms of contributions by the slow (s-) and rapid (r-) neutron-capture processes to the 
isotopic mix. While the result of the ID analysis of the Ba II feature indicates a 64±36% contribution 
of the r-process to the isotopic mix, the 3D analysis points toward a mere 15±34% contribution from 
this process, that is consistent with a solar-like mixture of barium isotopes. 

Keywords: hydrodynamics — line: profiles — stars: atmospheres — stars: abundances — stars: 
individual: HD140283. 



1 Introduction 

Elements beyond the Fe peak are predominantly pro- 
duced through successive neutron-capture reactions in 
two processes known as the slow (s-) and the rapid 
(r-) process. The distinction between the s- and r- 
process depends on whether the time-scale for neu- 
tron captures is longer or shorter, respectively, than 
the time-scale of radioactive decay of freshly synthe- 
sized unstable nuclei. Observational evidence and the- 
oretical studies have identified the s-process site in 
low- to intermediate-mass (~1.3-8 M©) stars in the 
asymptotic giant branch (AGB). In particular, the He- 
shell of thermally pulsating low-mass AGB stars is be- 
lieved to be the site of the so-called main s-process, 
which synthesizes nuclides heavier than Sr (see, e.g., 
IBusso et al.l Il999l ) . The r-process instead is usually 
associated with the explosive environment of Type II 
supernovae (SNell), although this astroph ysical site 
hasn 't been fully confirmed yet (see, e.g.. iSneden et all 
120081 '). 

Massive stars that end their evolution as SNell 
are shorter-lived compared with low-mass stars that 
eventually evolve into the AGB phase. One therefore 
expects the r-process to dominate the enrichment of 
the interstellar medium (ISM) in terms of heavy el- 
ements during the early stages of Galactic chemical 
evolutiorQ. According to this picture, the oldest most 
metal-poor stars should then only contain heavy ele- 
ments in the rela tive p roportions determined by the 
r-proccss (Truran 1981). 

Magain (1995) used a novel approach to infer the 
relative contributions of the r- and s-process spec- 
troscopically by looking at the detailed shape of the 
Ba II 4554 A resonance line. The odd Ba isotopes 



contribute to broaden the line and alter the symme- 
try of its profile via hyperfine splitting. As the r- 
and s-process produce rather different mixtures of odd 
and even Ba isotopes, the actual width and shape 
of the Ba ll 4554 A line are therefore dependent on 
the relative contribution of the two neutron-capture 
processes. Based oii the analysis of solar abundances 
(| Anders fc Grevessd IT989I ) bv lArlandini et al.1 l|l999t ). 
one expects a value of /o^d = 0.11 for the fractional 
abundance of the odd isotopes of bariurrQ, in case of 
a pure s-process mixture, and f^^^ = 0.46 in case of 
a pure r-process one. Magain challenged the Galactic 
chemical evolution scenario depicted above by deriv- 
ing a fractional abundance of odd Ba isotopes /odd = 
0.06±0.06 for the metal-poor subgiant HD140283 ([Fe/H] = 
—2.5): this value is fully compatible with a pure s- 
process production and excludes any significant con- 
tribution from the r-process. Magain's result cannot 
indeed be reconciled with the predictio ns of standard 
Galactic chemical evolution models (e.g. jTravaglio et al.l 
119991 ;) which indicate that the s-process signature in 

stars should become manifest only at higher metallici- 

ties ([ Fe/H]> — 1.5). More recentlv. iLambert fc Alle nde Prietol 
(|2002i ) re-observed and re-analysed the Ba ll 4554 A 
line in the spectrum of HD 140283 and derived, in con- 
trast to Magain, a fractional abundance of odd Ba 
isotopes /odd = 0.30±0.21, consistent with a pure r- 
proces s isotopic mixture. iLambert fc Allende Prietol 
(|2002l ') used a spectrum of superior quality than Mag- 
ain's both in terms of resolution (A/AA~200 000) and 
signal-to- noise ratio (S/N~550); their analysis, how- 
ever, was still based on ID model stellar atmospheres. 

Classical spectroscopic analyses based on ID model 
atmospheres, however, rely on a number of tunable 
fudge parameters. In particular, they cannot account 



^Assuming, of course, that SNell indeed are the site of 
the r-process 
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for Doppler broadening due to photospheric convective 
motions without introducing the micro- and macro- 
turbulence parameters. Moreover, convective flows also 
induce asymmetries and overall wavelength shifts in 
spectral line profiles, which cannot be reproduced by 
ordinary ID analyses. 

Three-dimensional hydrodynamical simulations of 
stellar surface convection, on the other hand, can self- 
consistently predict photospheric velocity fields and 
correlated t emperature and density inhomogeneities 
(jStein fc Nor dlund 1998; Asplund ct al. 1999); line shapes, 
asymmetries, and wavelength shifts can be accurately 
reproduced by using such simulations as 3D hydrody- 
namical model at mospheres, without relying on ad hoc 
free parameters (|Asplund et al. I I2OOOD . In the present 
contribution, we re-derive the fractional abundance of 
the odd Ba isotopes in HD140283 by means of both 
a ID and a 3D analysis o f the Ba II 4554 A line in 
the sp ectrum obtained by iLambert fc Allende Prietol 
(Hooa) at the W. J. McDonald Observatory (Mt Locke, 
Texas). In particular, we use a 3D model atmosphere 
of the metal-poor subgiant star to adequately disentan- 
gle the contributions of photospheric convective flows 
and hyperfine splitting to the broadening and asym- 
metry of the Ba 11 4554 A line. 



Table 1: Isotopic and hyperfine splitting compo- 
nents of the Ba 11 line at 4554 A. The compo- 
ne nts have been comp uted using the atomic data 



nents nave been comp uted using tne atomi c data 
bv lWendt et all (|198A and|vi Uemoes et al. 
for the splitting of the lower and upper energy lev- 
els. The oscillator strengths relative to the ^^^Ba 
component are given in the third column. The 
adopted loggf value for the line is +0.17. 



Isotope 


Wavelength (A) 


Relative strength 


^^4Ba: 


4554.0314 


1.0000 


"^Ba: 


4554.0003 


0.1563 




4554.0015 


0.1562 




4554.0019 


0.0625 




4554.0473 


0.4375 




4554.0500 


0.1563 




4554.0512 


0.0313 


i36Ba: 


4554.0317 


1.0000 


i37Ba: 


4553.9975 


0.1563 




4553.9986 


0.1562 




4553.9988 


0.0625 




4554.0498 


0.4375 




4554.0531 


0.1563 




4554.0542 


0.0313 


i38Ba: 


4554.0330 


1.0000 



2 Methods 

We synthesize the Ba 11 4554 A line profile with the 
help of both a ID and a 3D model atmosphere of 
HD140283, varying the Ba isotopic mixture to produce 
the best fit to the observed feature. We account for 
the isotopic and hyperfine splitting components of the 
Ba II 4554 A line listed in table [T] Isotopic wavelength 
shifts are actually too small to be resolved even with 
the highest-resolution spectrographs currently avail- 
able: in the present analysis, therefore, we efi'ectively 
do not make any distinction between the even isotopes 
as well as between the three odd isotopes. 

For the ID analysis, we employ a pla ne-parallel hy- 
drostatic LTE M ARCS model atmosphere (jGustafsson et aP 
I1975I : lAsplund ct al. 1997) of the subgiant star with 
the following stellar parameters: Tofi = 5690 K, \ogg — 
3.67 (cgs), and [Fe/H]= -2.50. Spectral line profiles 
are computed under the assumption of LTE. We stress 
that, contrary to the overall barium abundance, iso- 
topic abundances derived from the Ba 11 4554 A line 
are expected to be insensitive to departures from LTE 
and to the actual choice of stellar parameters. In 
the ID calculations, we adopt a micro-turbulence ^ = 

1.49 kms~^, based on the LTE a nalysis of the 51 Fe I 

and 13 Fe 11 lines given in Tab. 2 of lLambert fc Allende Prietol 
(|2002l '). With the above value, the iron abundance 
becomes independent from the equivalent width and 
chemical equilibrium is fulfilled (loge(Fe I) = 4.98±0.13 
and loge(Fe II) = 4.98±0.11). 

Before proceeding with the synthesis of the Ba 11 4554 A 
line we need to estimate the broadening from mecha- 
nisms other than hyperfine splitting. In order to do 
thi s, we fit the Fe lines from the s ame sample used 
by ILambert fc Allend e Prictc] lj2002t ) (Tab. 3 in their 
article) by assuming a rotational velocity Wrotsini — 
0.5 kms"'^ and varying the iron abundance, the cen- 
tral wavelength, and the FWHM of a Gaussian which 
we use toconvolve the line profiles. The convolution 
with a Gaussian accounts for macro-turbulent and in- 
strumental broadening (and, possibly, residual rota- 
tional broadening). We account for natural and lin- 
ear Stark broadening and we model coUisional broad- 
ening with neutral hydrogen atoms using the quan - 
tum mechanical calculations bv lBarklem et aP lj2000f ). 
From the analysis of the Fe lines, we derive an average 
value of 4.87±0.11 kms"^ for the FWHM of the Gaus- 
sian broadening. Using this estimate of the Gaussian 
broadening, we compute synthetic flux proflles of the 
Ba II 4554 A line for different barium isotopic mix- 
tures. We quantify the comparison between the theo- 
retical and observed profiles by means of a y'^-analysis 
similar to the one carried out bv I Asplund et al"] (jl999l ) 
to investigate ''Li/^Li ratios in metal-poor halo stars. 
We compute the according to the expression ~ 
^(Oi — Sif /a^ , with Oi and Si denoting the observed 
and synthetic flux at wavelength point i, respectively, 
and a is the inverse signal-to-noise ratio. The most 
likely isotopic mix is the one that minimizes x^ . 

For the 3D analysis, we use a simul a tion o f stel- 
lar surface convection by lAsplund et al. con- 
structed for the same stellar parameters as for the ID 
CcLS6, £LS 3, time-dependent 3D hydrodynamical model 
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Figure 1: Upper panel: Synthetic {eontinuous 
black line) versus observed (grey symbols) profiles 
of the Ba ii hue at 4554.0 A for the best fitting 
Ba isotopic mix in the ID analysis. Lower panel: 
Relative difference between observed and synthetic 
profiles. 



atmosphere. We select a 35-minute long sequence of 
30 simulation snapshots sampled at regular time in- 
tervals and downgrade the resolution of the original 
simulation to 50x50x82 grid-points prior to the spec- 
tral line formation calculations. We compute LTE flux 
profiles in 3D by solving the radiative transfer equa- 
tion along 33 rays (four /x-angles, eight 0-angles, and 
the vertical) for all grid-points at the surface, perform- 
ing then a disk integration and a time average over 
all snapshots. As we already account for the effect 
of Doppler shifts induced by the velocity fields in the 
3D simulation, we ignore altogether the micro- and 
macro-turbulence parameters. We then determine the 
amount of broadening contributed by processes other 
than hyperfine splitting: we fit the profiles of the Fe 
lines similarly as in the ID analysis above but varying 
the rotational velocity instead of the Gaussian broad- 
ening, which, in 3D, we keep fixed and equal to the 
instrumental broadening (corresponding roughly to a 
FWHM= f.5 kms~^, as estimated from the resolu- 
tion of the observed spectrum). Natural and coUisional 
broadening are included and treated in the same way 
as in the ID calculations. From the 3D analysis of 
the Fe lines, we estimate an average rotational veloc- 
ity of 2.58±0.30 kms"\ We then follow the same x^- 
minimization procedure as in the f D analysis to deter- 
mine the most probable barium isotopic mix from the 
fitting of the observed Ba ll 4554 A line profile. 



3 Results 
3.1 ID analysis 

Figure [T] shows the best-fitting synthetic profile to the 
observed Ba ll 4554.0 A line in the ID case. The de- 
rived fractional abundance of odd Ba isotopes in fD 
is /odd = 0.33, which suggests a contribution predom- 
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Figure 2: Thick continuous line: Resulting Ax^ — 
- Xmin fro™ the fit of the Ba ii 4554 A line 
as a function of the fractional abundance of odd 
Ba isotopes for the combined Gaussian (FWHM= 
4.87 kms"^) of the mean macro-turbulence and 
instrumental broadening in the ID case; also 
shown are the cases for FWHMiAFWHM {thick 
dashed and dot- dashed lines). A value Wiotsini = 
0.5 kms~^ is adopted for the rotational velocity. 
Thick grey lines: fractional abundances of odd Ba 
isotopes in case of pure s-process and r-process 
isotopic mixes (S. Bisterzo 2008, private commu- 
nication). 



inantly from the r- 
ium isotopic mix. 



process (namely 64%) to the bar- 
Figure [5] shows the Ax^ = ~ 
Xmin curves considered in the x^-minimization pro- 
cedure. The la, 2a, and 3a confidence limits cor- 
respond to Ax^ = f, 4, and 9, respectively (indi- 
cated by the dotted lines in the figure). The half- 
width of the /odd interval for which Ax'^ < f is a 
measure of the uncertainty in fractional abundance 
of odd isotopes due to the finite signal-to-noise; for 
S/N= 550 the uncertainty on /odd is ±0.023. Fig- 
ure [2] also shows how the position of the x^ mini- 
mum depends on the Gaussian broadening: a change 
of zhO.lf kms~^ translates into a change of =f0.06 in 
/odd- The fractional abundance of odd Ba isotopes is 
also sensitive to the choice of micro-turbulence param- 
eter (which affects the Gaussian broadening determi- 
nation) and to the adopted surface gravity (which af- 
fects pressure broadening of the lines). For the present 
analysis, we adopt the estimates o f =f0.11 and =f0.01 by 
iLambert fc AUende Prietd (|2002t ) of the errors on /odd 
due to the uncertainties on micro-turbulence (±0.2 kms~ 
and \ogg (±0.2 dex), respectively. Our estimated total 
uncertainty on /odd is then A/odd ~ 
VO.O232 ± 0.062 _^ q \i2 0.012 -0.13 which directly 
corresponds to an uncertainty in the contribution of 
r-process fraction of Ba isotopes of about ±36%. 
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Figure 3: Same as Fig. [3] but for the best fitting 
Ba isotopic mix in the 3D analysis. 



3.2 3D analysis 

The best-fitting synthetic Ba ll 4554.0 A line profile in 
the 3D analysis is shown in Fig. |3] It is clear from the 
comparison of the residuals, that the synthetic profile 
in 3D gives a significantly better fit to the observed 
spectral line than the synthetic ID profile. The esti- 
mated fractional abundance of odd Ba isotopes in the 
3D analysis is /odd = 0.15, which — contrary to the ID 
case — points toward a predominantly s-process iso- 
topic mix (only 12% contribution from the r-process). 
Figure U] shows the Ax^ ~ ~ Xmin curves con- 
sidered for the x'^-minimization procedure in the 3D 
analysis. Similarly as for the ID case, we use the in- 
formation in the figure to quantify the uncertainties 
on /odd due to the finite signal-to-noise and the error 
on the estimated rotational velocity: these are 0.024 
and 0.12, respectively. For the error on /odd due to 
the uncertainty on surface gravity, we adopt the same 
value as in the ID case. We caution that a change 
in surface gravity in practice also affects the velocities 
predicted by the simulation. In this preliminary 3D 
analysis of the Ba ll 4554.0 A line, however, we ne- 
glect this aspect; we defer the study of the effects of 
changes in stellar parameters on the estimate of the 
uncertainty on /odd to a future and more comprehen- 
sive work. Finally, we do not consider any contribu- 
tion to the error on /odd from micro-turbulence, since 
this parameter is absent in the 3D analysis. The es- 
timated total uncertainty on /odd in the 3D analysis 
is then A/odd = VO.0242 0.122 -h 0.01^ ;~.o.l2 which 
directly corresponds to an uncertainty in the contribu- 
tion of r-process fraction of about ±34%. 



4 Discussion 

The fractional abundance /odd = 0.33±0.13 of odd Ba 
isotopes derived in the present ID analysis is in excel- 
lent agreement with the r esult obtained by 
iLambert fc Allende Prietol l|2002D . Contrary to the find- 
ing of lMagainl(|l995l l. our ID estimate of /odd indicates 
a large contribution from the r-process (64%) to the 
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Figure 4: Thick continuous line: Resulting Ax^ 
from the fit of the Ba ii 4554 A line as a func- 
tion of the fractional abundance of odd Ba isotopes 
for a rotational velocity Wrotsini = 2.58 kms~^ 
in the 3D case; also shown are the cases for 
t^rot sin i± At;rot sin j {thick dashed and dot- dashed 
lines). A value FWHM= 1.5 kms^^ is assumed 
for the instrumental broadening. 



barium isotopic mix in HD 140283 and seems to exclude 
at a 3(T level a pure s-process contribution. 

Interestingly, the 3D analysis of the same feature 
leads to a radically different result, namely that the 
Ba II 4554.0 A line profile is best-fitted assuming in- 
stead a solar-like barium isotopic mix with only a 15% 
(±34%) contribution from the r-process. Moreover, in 
3D, a pure r-process isotopic mixture can be excluded 
at a 2a level from the analysis. The reason for the 
difference between the results of the ID and 3D analy- 
ses is essentially due to the differences between the ID 
and 3D modelling of line broadening and asymmetries 
as sketched in Sec. [2] For instance, in the 3D calcula- 
tions, the simulation's velocity fields naturally induce 
a "C-shaped" asymmetry in the flux profile of the line. 
This implies that a different fraction of odd barium iso- 
topes is necessary in 3D comparing to the ID case to 
properly model the shape of the Ba ll 4554.0 A line 
profile. We would like to caution at this point the 
reader that the our result for HD 140283 cannot be ap- 
plied straightforwardly to the ID analysis of the Ba ll 
4554.0 A feature in the spectra of other stars. The 
sign and magnitude of the 3D— ID correction to the 
derived fractional abundance of odd isotopes may de- 
pend in general on the strength and detailed shape of 
the observed Ba ll 4554.0 A profile. 

Our 3D analysis of the Ba ll 4554.0 A feature in 
the spectrum of HD 140283 seems unable t o settle down 
the controversy raised by iMagainI ()l995l l in favour of 
standard Galactic chemical evolution scenarios. On 
the contrary, our result seems to corroborate the pos- 
sibility that the metal-poor star HD 140283 possesses 
a strong s-process signature in terms of isotopic Ba 
abundance against the expectations from theoretical 
models of Galactic chemical evolution. We would like 
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to stress however that, although our and Magain's re- 
sult are apparently in agreement with each other, our 
finding is based on rather different premises. First, we 
rely on a higher-resolution and higher signal-to-noise 
spectrum. Second, we use a 3D hydrodynamical model 
atmosphere of the star to synthesize the profiles of the 
Ba II 4554.0 A feature and other Fc lines, which im- 
plies that our modelling of non-thermal broadening is 
more robust and, contrary to the ID analysis, inde- 
pendent on the tunable micro- and macro-turbulence 
parameters. 

As a final note, we cannot completely rule out 
the possibility that HD140283 has been polluted by 
s-proccss material from an AGB star and may not be 
representative of halo stars at the same metallicity in 
terms of isotopic Ba abundance. Further investigation 
is therefore necessary to draw any significant compar- 
ison with Galactic chemical evolution models. We in- 
tend to extend in an upcoming paper the 3D — ID anal- 
ysis of the Ba ll 4554.0 A feature to a larger sample of 
halo stars and study in particular the barium isotopic 
abundance as a function of metallicity. 
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